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PREFACE 


The report contains the results of several papers 
related to modeling using a class of the bivariate gamma 
distribution. The separate papers contain loosely related 
subjects pertaining to this problem. Since the separate 
papers were prepared at different times during the contract 
period and have been submitted for publication in the open 
literature and each paper is intended to be self-contained, 
there is some redundancy in tables and illustrations. 

Each of the papers in this report were extensions 
and/or generalizations of the results given in NASA TM- 
82483, entitled "A Bivariate Gansna Probability Distribution 
with Application to Gust Modeling," by 0. E. Smith, S. I. 
Adelfang, and J. D. Tubbs. A modification of this paper is 
currently under review by Communications in Statistics ■ 

The first paper in this report, entitled "A Note on 
the Ratio of Positively Correlated Gamma Variates," has 
been accepted for publication in Communications in Statis- 
tics and it presents some new analytical results using a 
class of the bivariate gamma distribution. Comparable 
results were available in the open literature using a 
different class of the bivariate gamma. 

The second paper is entitled "A Method for Determin- 
ing if Unequal Shape Parameters are Necessary in a Bivariate 
Gamma Distribution" and is an application of the results 
given in the first paper and addresses questions concerning 

ii 


hypothesis tests for equality of shape parameters from 
correlated gamma distributed variates. This paper is 
currently under review by Technometrics . 

The third paper, entitled "A Differential Equations 
Approach to the Modal Location for a Family of Bivariate 
Gamma Distribution," contains extensive analytical results 
for the location of the mode as a function of the free para 
meters. To the authors' knowledge this is the only such 
representation for a non-gaussian bivariate distribution. 
This paper has been submitted to SIAM on Scientific 
and Statistical Computing . 

The fourth paper is a report summarizing the analysis 
of some wind gust data using the analytical results devel- 
oped in relationship to the modeling application. 
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A NOTE ON THE RATIO OF POSITIVELY CORRELATED GAMMA VARIATES 


J. D. Tubbs 

Department of Mathematical Sciences 
University of Arkansas 
Fayetteville, Arkansas 


0. E. Smith 


Systems Dynamics Laboratory 
NASA Marshall Space Flight Center 
Huntsville, Alabama 


ABSTRACT 

Mielke and Flueck (1976) derived the density function and 
corresponding moments for the ratio of correlated gamma distrib- 
uted variates. They considered a class of bivariate gamma dis- 
tributions suggested by Cherian (1941) and David and Fix (1961). 
Recently, Lee, Holland, and Flueck (1979) derived some additional 
distributional results using this class of functions. This paper 
derives similar results using a different class of bivariate 
gamma distributions. 


1. INTRODUCTION 

Mielke and Flueck (1976) derived the distributional results 
for the ratio, R, of correlated gamma distributed variables. 

There are several classes of the bivariate gamma distribution [three 
are summarized in Mardia (1970) and an additional two in Johnson and 
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Kotz (1972)]. Mielke and Flueck (1976) derived the distributional 
results for the ratio, R, of correlated gamma distributed variables 
using the Cherian-David-Fix class of bivariate gamma random variables 
[Cherian (1941) and David and Fix (1961)]. That is, let X, Y, and P 
denote independent gamma random variables with common scale parameter 
X and respective shape parameters a -C, 6 “ S» and £, for 0 < £ < 
min(a ,6) . 

Then it can be shown that the bivariate probability density 
function for U X + P and V ■ Y +.P is given by 

r min (u,v) 


f u,v (u » v) 


o pQ ~ 




( 1 . 1 ) 


for K - r(a-5)r(6-5)r(C) 

when the scale parameter X is assumed to be unity. Mielke an' 
Flueck (1976) showed that (1.1) can be written as 


f U,V (u ’ v) 


a-1 B-S-l -(u+v) 

, u v e 

r(a)r(3-0 


, a-C-1 B-l -(u-W) 
1 u v e 


F^Cfl+C-usBjv/v.-u) if 0<u<v 


( 1 . 2 ) 


rCa-Or(s) 


F£(S,H£-a:8,v/u,-v) if 0<v<u 


m n 

where F*(a,b,c:x,y) =■ £ -t-z j^y^y x is a "degenerate" 

m,n a 0 ' 'n+n 

two variable hypergcomctric function [Gradshtcyn and Ryzhik 
(1967), p. 1067] and (a) n “ F(a+n)/r(a). Thus, U and V are gamma 
random variables with shape parameters a and B and positive depend- 
ence parameter £. In particular, E(U) “ Var(U) = a, E(V) * Var(V) 

- B, and Cov(U,V) •> £. 

Mielke and Flueck (1976) derived the density function for 
R ** U/V using a change of variables. That is, 
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r° ^l+r)^" 01 ^ F^.a+B-C.l+S-B.atr/l+r.r) 


I B(o-£,3) 


f R (r) - { 


if 0 < r < 1 


(1.3) 


r ^* -1 (l+r) C ~ a ~ g F.(C.ct+B-C. l+?-a,B :l/l+r,l/r) 
B (a-C,B) lf X< r 


• (a)_^(b) (c) 

where F 1 (a,b,c,d:x,y) - I tt ” n i x m y n , |x| < 1 |y| < 1 

m,n=0 m+n 

is a two variable hypergeometric function [Gradshteyn and Ryzhik, 
(1967), p. 1053], and B(a,b) - r(a)r(b)/r(a+b). In addition, they 
show that the integral moments of R are given by 

„ S fo-O ,(€)_, 

E(R S ) - jEq (p (S-j) " 1 for s >, 0. (1.4) 

s 


In particular, 

E(R) ° «B-C/6(B-1), 3 > 1 


(1.5) 


EfR 2 ! o + 2 (a-e)e . (a-g+1) (a-g) . , 

K } 3(3+1) 3(3-1) (B-1HB-2) » p i 


(3-1) (3-2) 


Recently, Lee, Holland, and Flueck (1979) were able to obtain 
comparable results for density of R using the Cherian-David-Fix 
class of densities by expressing f 0 as a weighted difference of 
hypergeometric functions. The purpose of this paper is to derive 
comparable results for R using a different class of the bivariate 
gamma distribution. This class is a special case of the one 
suggested by Jensen (1970) as modified by Gunst and Webster (1973). 
The next section contains a brief discussion of this class of 
distributions. In section 3 the derivation of f is given using 
this class of functions. Section 4 outlines a possible applica- 
tion for the probability function in the area of hypothesis 
testing for the equality of shape parameters in the presence of 
correlation. 
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2 . GUNST AND WEBER CLASS OF BIVARIATE GAMMAS 


Gunst and Weber (1973) proposed a computationally feasible 
method for deriving the joint density function for the bivariate 
chi-square distribution. Since the chi-square is a special case 
of the gamma, this method was used for the bivariate gamma case. 
That is, a bivariate gamma density function for U and V with 
common scale parameter X * 1 and shape parameters a, 3, (a < 3) 
is given by 


f(u,v) 


u a-l v 3-l e -[ (u+v)/(l-n)] " 
(i-n) a r(a)r(B-a) j 


£ n~* +k T (3-a+k) (uv)^v k 
k (l-nJ^^CS+j+kJjlkl 


where n “ p /(3/a) , p is the correlation coefficient between the 
variables U and V. Gunst and Webster (1973) suggested this class 
of densities in that they are computationally tractable and do 
not involve mathematical functions, such as Laguerre polynomials 
or convoluted sums [Jensen (1970) and Kibble (1941)]. Smith, 
Adelfang, and Tubbs (1982) discuss this class of densities in 
greater detail. 

In the next section the distributional properties for the 
ratio, R, are derived using the Gunst-Webster class of bivariate 
gammas. 

3. RATIO OF CORRELATED GAMMA VARIATES 

By letting R = U/V and S = U+V, the joint pdf for R and S 
can easily be shown to be 


f R,S <r ’ S > 


£ £ 

i * 


(1+r) 2 



a+j-1 



S+j+k-1 

e 


-s/(i-n) 


(3.1) 


where Cj •= [ (l-n) a r (a) (B-a)] -1 , c. 


n j+k r(s-a+k) 


2 (l-n) 2j+k T(3+j+k)jlk! 

by integrating over S the pdf for R becomes 


Hence, 
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QO CO 

f K (r) - (1-1,)®^ CjC t B( n *j, 6 +J+k)" 1 r‘‘ +J ' 1 /(l+t) a4<+2 -!' fk (3.2) 

" her ° ~-l ’ « k " (6-c) k n k /k!. <a) n - «,«)«.), «„ d 

B(a,b) - r(a)r(b)/r(a+b). 

Whenever the shape parameters are equal then the density func- 
tion for R is given by 

f R (r) - (1-n) 0 ^ Cj B(« + j,a + j)- 1 r 0+ > 1 /(H r) 2^ (3#3) 

From (3.2) it can be shown that the m th raw moment for R is 
given by 

E(R ) - (1-n) c k B(a+j+m,B+j+k-m)/B(a+j,B+j+k) (3.4) 

if m < B. In which case, it follows that 

CD CO 

E(R) - ( 1 -n) e E Cj E c t (a+j)/ ( 8 +j+k-l) (3 . 5) 

j 1c 

and 


E(R ) - (1-n) 0 Z Cj Z c k (a+j) (a+j+1)/ (B+j+k-1) (B+j+k-2) (3.6) 

j k 

Whenever n ° 0, then 

E(R) - a/ (B-l) , E(R 2 ) = a (a+1)/ (B-l) (B-2) (3>7) 

which agrees with the values given by Mielke and Flueck (1976) 

whenever E - 0 and wxth Lee, Holland, and Flueck (1979) whenever 
a “ 0. 

Lee, Holland, and Flueck (1979) discuss some of the mathe- 
matical properties for the density of R for various values of 
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a,0, and n. They demonstrated that the density can be » at r-1 
whenever either of the shape parameters is less than one. How- 
ever, in the Gunst-Webster construction by assuming that o > 1 
and a < 0 the density function given in equation (3.2) is stable. 
Figures 1-4 illustrate the various shapes that f R (r) has as a 
function of the three parameters. 

A definite computational advantage of equation (3.2) versus 
equation (1.3) stems from the ability to compute the tail proba- 
bilities for R. By letting a=a+j and b=0+j+k, we have 


W ' «-»> 9 * ‘J l 'k ='tt 2a , 2b i b V aJ 


(3.8) 


where F.^ denotes a random variable from an F-distribution with 
r and s degrees of freedom. Note if n =0, then (3.7) becomes 


W = PfF 2n,20 i e V Q] 


(3.9) 


which agrees with the well known results concerning the ratio of 
independent chi-squares. Furthermore, if n t 0 and a = 0 then 
(3.7) becomes 


F R (r 0 } 


(1-n) I c p[f. 

i j 


2 (a+j),2(a+j) = r Q^ 


(3.10) 


j which is similar to an expression given by Johnson and Kotz (1972), 

- Chapter 40, Section 3. 

t 4. APPLICATION 

J 

In this section an application is given for computing the 

cdf of R, given by equation (3.7). Diagram 1 defines the area 
given in equation (4.1). 
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By letting cot 0 Q - U/V - r Q and C(0) - l-F R ( r(J ), one haa 
G(e) - (l-n) e Z Cj Z c k P[r 2b ^ 2a < (a/b)tan 0] (4.1) 

Figure 5 contains the graph of the function G(6) versus 0 for o * 1 
and B - 1, 2, or 3 and n - 0, .25, .50, and .75. From this figure 
and other cases which are not included one observes that whenever 
a ■ 0 then C(45°) - .5 and G(45°) < .5 whenever a < 0. This obser- 
vation and additional properties were used in developing a test for 
the hypothesis 


H Q : o - 0 vs. H a : a < 0 


(4.2) 


The procedure is presented in Tubbs (1983) and uses the Craaer- 
Von Miseo criteria for testing (4.2). That is, define 

\ - " / l V r > - F,«> 2dF R M (4.3) 

where F R (r) is the cdf for the null distribution given in (3.10). 

F n (r) is the empirical distribution for ** u^/v^ an< * t * ie r^'s 

are arranged in increasing order. Whenever H n is true, then YJ is 

u n 

distribution free and has a convenient computational fora given by 


W 


1 

12n 


+ Z IZ - 
i-1 1 2n 


(4.4) 


where * F R (r^). ^0 rejected if exceeds a specified 

critical point. Tubbs (1983) considers the properties of this 
test procedure in greater detail. 


5. CONCLUSIONS AND SUMMARY 

This paper derives both the density and the distribution func- 
tions for the ratio of positively correlated gamma variates using 
a codification of Jensen's bivariate gamma distribution. The 
expression for the moments differ from those given by either 
Miclke and Flueck (1976) or Lee, Holland, and Flueck (1979). 
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However, all the expressions are Identical whenever the variates are 
uncorrelated. A principal advantage found in this representation 
stems from the ability to compute the CDF of the ratio. The value of 
the CDF for the ratio was shown to have potential application to the 
problem of testing for equality of shape parameters in a particular 
family of the bivariate gamma distribution. 
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A METHOD FOR DETERMINING IF UNEQUAL SHAPE 
PARAMETERS ARE NECESSARY IN A BIVARIATE 
GAMMA DISTRIBUTION 

J. D. TUBBS 

Department of Mathematics 
University of Arkansas 
Fayetteville, Arkansas 

ABSTRACT 

A procedure for aiding an experimentalist in decid- 
ing between four and five parameters in a Jensen's type 
bivariate gamma distribution is presented. The procedure 
is based upon the properties of the CDF for the ratio of 
correlated gamma distributed variates. The criteria of 
interest is posed in a test of hypothesis setting and 
results are presented using the Cramer-Von Mises test 
of fit. 

1. INTRODUCTION 

Smith and Adelfang (1981) discuss the applicability 
of a bivariate gamma distribution as a parametric model 
for wind gust amplitude and length. In modeling this 
bivariate data with a gamma distribution, it was neces- 
sary to find a distribution that would allow for correla- 
tion between the random variables X and Y when the marginal 
distributions are univariate gammas with possibly unequal 
shape and scale parameters. That is, X - G(y v ,B v ) and 

A* X 


Y “ GCy^.Py) where the probability density function for 
Z - G (y , B ) is given by 

f z (z) » 0 v z Y " 1 e“ 0z /r ( y). (1.1) 

A brief survey of the open literature reveals that 
there are several classes of the bivariate gamma distribu- 
tion. One need only consult Mardia (1970) and Johnson and 
Kotz (1972) to find five classes of the bivariate gamma 
distribution [Kibble (1941), Cherion (1941), McKay (1934), 
Jensen (1970), and Moran (1969)). Of these classes only 
Jensen (1970) and Moran (1969) allow for unequal shape 
parameters and both of these have computational limitations 
which affect their utility to the experimentalist. Recently, 
McAllister, nee, and Holland (1981) and McAllister (1933) 
have addressed the limitations with Jensen's model and 
provided results which overcome many of the computational 
difficulties. However, at the time of Smith et al. (1983) 
development these results were not available. Hence, they 
modified a bivariate chi-square model given by Gunst and 
Webster (1973). This allows for possibly unequal shape 
parameters and is computationally tractable. The model is 
not as general as that given by Jensen (1970), however, one 
can derive the bivariate model given by Kibble (1941) as 
a special case whenever the shape parameters are equal. In 
this paper, the unequal shape parameter model will be 
referred to as the five-parameter model and the equal shape 
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case as the four-parameter model. Smith, Adelfang, and 
Tubbs (1983) discuss the properties of these distributions 
and it is apparent that the four-parameter has numerous 
computational advantages over the five-parameter model. 

So if one assumes that the data is correctly modeled by 
this class of the bivariate gamma distribution, a question 
of practical interest becomes, How does one decide if the 
five-parameters are really necessary? The purpose of this 
paper is to present a procedure which would aid the exper- 
imentalist in answering the above question. The problem 
is posed in a hypothesis testing setting. That is, test 
the hypothesis 


H n = Y, 


versus 


( 1 . 2 ) 



'x 



(1.3) 


It should be noted that the proposed method is not an 
omnibus test of fit for the bivariate gamma against all 
other possible models. Instead the procedure is intended 
for deciding between the four or five parameter models as 
given in Smith, Adelfang, and Tubbs (1983). 

The next section contains the distributional results 
needed for the test of hypothesis (1.2). The test proce- 
dure is given in section 3 and evaluated in section 4. 
Section 5 contains a summary and remarks concerning some 
of the limitations of the procedure. 
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2. DISTRIBUTIONAL RESULTS 
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Smith, Adelfang, and Tubbs (1983) modified a bivariate- 
Chi square distribution given by Gunst and Webster (1973) 
and obtained the density function given by 

f(x,y) = K,/K 2 2 i c..x j (ny) j+k (2.1) 

1 z j=0 k=0 JK 

where 

= x x y Yy exp-{ (x+y)/(l-n)} , 

K 2 = (l-n) Yx r(Y x )r( Yy - Yx ), 

c jk - n j+k r(Y y -Y x +k)/{(l-n) 2 j+k r( Yy +j+k) jl kl), 

and x = X6 , y - Yf , B , B„ are known scale parameters, 
x y x y 

n “ p Y y /Y x » P is the correlation coefficient between the 
variables X and Y. The joint probability distribution func- 
tion is given by 

F(x , y ) = Pr[X < x , Y < y ] 

= J £ £ d., H(y +j ,x /(1-n)) 

j=0 k=0 Jk x o 

•H( Yy +j+k, y Q /(l-n)) (2.2) 

where 

J = (1-n) y /r (y x )t (Y y -Y x ) . 

d jk ” n ^ +kr (Yy"Y x +k)/r (Y y +j+k) jl k! 

H(a,x) = /“t a_1 e _t dt. 



Equations (2.1) and (2.2) are for the unequal shape para 

meters and will be referred to as the five-parameter model. 

It should be re-emphasized that this model is not completely 

general in that one assumes that y„ > Y v and the correlation 

y x 

between variables X and Y are restricted to the interval 
[ 0 , n^Yjj/Yyl for ne[0,l]. 

If Y v “ Y v " Y then it can be shown that (2.1) and 

x y 

(2.2) reduce to the well known functions given by Kibble 

(1941). • That is, the density function is given by 

f(x.y) = (xy) Y ' 1 exp-{(x+y)/(l-n)>/r(Y) 

• £ (nxy)/(l-n)' ! ) j /r( Y +j) jl (2.3) 

j=0 

and the distribution function becomes 

F(x,y) = (l-n) Y /r(Y) i n^/r(y+j) j 1 

j=0 

*H(y+] ,x/ (1-n)) H(Y+j,y/(l-n)). (2.4) 

Equations (2.3) and (2.4) will be referred to as 
the four parameter model. A comparison of the distribution 
function given in (2.2) and (2.4) reveals that there are 
distinct differences in terms of the computational com- 
plexity. Thus for computational reasons the experimental- 
ist would like to know how much greater does Y y have to 
exceed y x before equation (2.2) is really necessary. 

Ideally he would like to answer this question before using 
both (2.2) and (2.4) then selecting the results which are 
more gratifying. In order to address this issue, this 
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paper considers the problem of testing hypothesis (1.2) 
versus (1.3) using an univariate random variate given 
by the ratio of X to Y, R - X/Y. Tubbs and Smith (1983) 
derive the density and distribution functions for R when- 
ever the bivariate density is either (2.1) or (2.3). That 
is, if equation (2.1) holds then the density function for 
R is given by 

f R ( r ) = (1-n) y ^“^CjC^B(a,b) -1 r a ~ 1 / (l+r) a+b ( 2 * 5 ) 

where B(a,b) = r (a)r (b)/r (a+b) , Cj - (a) jn j /jl, 

C k = (b-a) k n k /kl, a = Y x +j , b = Y y +j+k, and (a) R - r(a+n)/r(a). 
The distribution function for R is given by 

F R (r o> = Pr f X / y 4 r 0 i 

* k I 0 C j C k Pr[F 2 a .2b i b *o <2 -6) 

where F^ s denotes a random variable from an F-distribution 
with r and s degrees of freedom. The correspo ding functions 
whenever y x ~ Y y = y are given by 

f R (r) = (l-n) Y i c. B(a,a) -1 r a " x /(l+r) 2a (2.7) 

j=0 J 

and 

F R (r 0 ) = (1-n)^^ Cj Pr[P 2a 2a < r 0 ) (2.8) 


where a = y + j . 
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3. HYPOTHESIS TESTING 


Since R = X/Y is a univariate random variable it is 
informative to graph Fj^(r) versus r. However, since r > 0 
a more meaningful graph can be produced by letting 9 =* 
cot~^r and G(e o ) - l-F r (r Q ) where 9 Q = cot“^r o . The area 
corresponding to F r (r Q ) is shown in diagram 1. Further- 
more, it follows that 

G(9 ) = (1-n) y Z 2 c.c k Pr[F ? 2h < (a/b) tane] 
j=0 k=0 J K za ,^ D ° 

(3. 1) 

in the’ five-parameter model and 


G(9 ) = (1-*0 Y 2 c Pr[F„ 
° j=Q J 


2a 


4 tan9 Q ] 


(3.2) 


in the four-parameter case. 

Since 9 is restricted to the finite interval (0, tt/2), 
it is somewhat instructive to plot G(e) versus 9 as func- 
tions of the free parameters, y x , y y and t\. As in Tubbs 
and Smith (1983) the scale parameters are assumed to be 
known and hence equal to one. This restriction will be 
addressed later in the paper. Figures 1-3 contain some of 
the illustrative cases. 

From these plots one observes that G(45°) ** .5 when- 
ever the four-parameter model holds and G(45°) < .5 in 
the five-parameter models. Rather than just using this 
observation a function was selected to measure the distance 
between these distribution functions. The Cramer-Von Mises 








6 angle in degrees 


FIGURE 1. F p (9) vs. 0 
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type goodness-of-fit procedure was selected since the tost 
is distribution free whenever the parameters are specified. 
Furthermore the test statistic is easy to compute. 

Let 

m/2 - 

W n " n / < G (°> ' G n (°)> dG(0) (3.3) 

O 

where G(e) is given in (3.2), G n (o) is the empirical dis- 
tribution function of tan"^(r^). r^ “ are 

arranged in increasing order. Whenever hypothesis (1.2) is 
true, then W n has the convenient computational form given 
by 

W„ - l/12n + s {z. - (3.4) 

n i=l 1 /n 

where = G(e^). Furthermore, from Anderson and Darling 
(1951) one can reject (1.2) whenever W n exceeds a specified 
critical point. These critical points are given from Ander- 
son and Darling's asymptotic distribution. Stephen (1976) 
defines a procedure for modifying the critical points for 
small samples, however, the underlying problem of modeling 
bivariate data will probably dictate large sample sizes. 

4. EVALUATION OF THE TEST PROCEDURE 

In this section the procedure defined in the previous 
section is evaluated. The evaluation is performed in two 
parts. The intent of the first part was to determine 
whether or not the procedure even works. That is, are the 
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apparent visual differences between the function G(o) as 
seen in Figures 1-3 significant in the "Cramer-Von Mises" 
metric. The second part of the evaluation concerns the 
robustness of the procedure to the nuisance parameters. 

In the first part, let 

tt/2 « 

d_< 6) - n / {G(o) - A(e)rdG(e) (4.1) 

o 

where G(e) is given in (3.2) and A(o) is given by (3.1) 

when Y v " Y v + <5, for & >0. For positive integers n, 

y x 

compute . 

a n («) “ p r [W n > D n (fi)i- (4.2) 

If the alternative hypothesis given by 

H l : < Yy = Y x + « (4.3) 

holds, then the expected value of W in (3.3) is given by 
D n (£). Hence, a n (6) is the expected type I error of test- 
ing hypothesis (1.2) as a function of 6. Table 1 contains 
the value of <* n (<5) for various values of the parameters. 
The a n (6)'s were computed using Tiku's approximation to 
the asymptotic distribution of W n [Tiku (1965)]. 

For example, from Table 1 one would expect the test 
to reject integer (6=1) differences between the shapes for 
X and Y at the 95% significant level whenever n > 50. 

The procedure used to generate the values in Table 1 
is somewhat unconventional, however, they do indicate that 


1 

1 

J 

i 

» 
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the test procedure would be sensitive to differences in 
the shape parameters that exceed unity. A Monte Carlo 
simulation was also performed. The results are not 
reported in the interest of space and since the simulation 
was quite limited. A detail simulation is very expensive 
due to the computational costs in computing the null dis- 
tribution G(e) needed in evaluating type I errors. It 
is especially costly to simulate any type II errors. In 
spite of these restrictions upon the simulation's merit, 
the results were supportive of the expected results given 
in Table 1. 

The second part of the evaluation is concerned with 

the question of robustness of the test to the unspecified 

parameters, namely, p and g , g . In order to determine 

x y 

the sensitivity of the test to the misspecified correla- 
tion coefficient p, the following distance was evaluated 

for different values of y = 

x y 

D n (p) » n 1 /2 {G(Q) - B(G)} 2 dG(e) (4.4) 

where G(o) is given in equation (3.2) when p = 0, and 
B(o ) is given by equation (3.2) whenever p > 0, for 
p = .25 (.25). 75. Table 4 contains the type 1 errors 
a n (p) given by 

p r[W n > D n (p ) J - a n (p) (4.5) 

for different values of n and y — y *= y . a_(p) in 

x y ii 
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Table 1 . Tail Probabilities for “ n CO * 


n _ 

n 

6^.25 

.50 

.75 

1.00 

1.25 

1.50 

1 . 75 _ 

2.00 

1 0 

20 

1.00 

.45 

.22 

.12 

.06 

.04 

.02 

.02 


50 

.53 

.12 

.02 

.01 






100 

.25 

.02 

.01 






.25 

20 

1.00 

.41 

.18 

.09 

.05 

.03 

.02 

.01 


50 

.48 

.09 

.02 

.01 






100 

.21 

.01 







.50 

20 

.87 

.33 

.13 

. 96 

. 03 

.02 

.01 

.01 


50 

.41 

.06 

.01 







100 

.15 

.01 







.75 

20 

.69 

.21 

.07 

.02 

.01 

.01 




50 

.27 

.02 








100 

.06 








2 0 

20 

1.00 

.82 

.48 

.27 

.17 

.13 

'.07 

.04 


50 

1.00 

.36 

.13 

.04 

.02 





100 

.60 

.13 

.02 

.01 

.01 




.25 

20 

1.00 

.72 

.40 

.22 

.13 

.07 

.04 

.03 ' 


50 

.91 

.29 

.09 

.01 






100 

.52 

.09 

.01 






.50 

20 

1.00 

.59 

.30 

.15 

.08 

.04 

.03 

.01 


50 

‘ .73 

.21 

.05 

.02 






100 

.41 

.05 

.01 






.75 

20 

1.00 

.40 

.16 

.06 

.03 

.01 




50 

.52 

.09 

.01 







100 

.24 

.02 







3 0 

20 

1.00 1 

.00 

.66 

.43 

.28 

.18 

.12 

.08 


50 

1.00 

.55 

.25 

.11 

.04 

.02 

.01 



100 

.87 

.26 

.07 

.01 





. 25 

20 

1.00 1 

.00 

56 

.35 

.21 

.13 

.08 

.05 


50 

1.00 

.46 

.18 

.07 

.02 

.01 




100 

.59 

.11 

.01 






.75 

20 

1.00 

.56 

.26 

.12 

.05 

.02 

.01 



50 

.72 

.18 

.04 

.01 






100 

.39 

.04 

.01 







*if n n (0 < . 01 , then the entry is left blank . 
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(4.5) is the expected type I error as a function of the 

nuisance parameter p. It should be mentioned that the 

distribution for W n is not the same as that given by 

Anderson and Darling asymptotic approximation since the 

nuisance parameter p is unspecified [cp. , Stephen (1976)], 

however, it does not appear feasible to follow Darling's 

procedure for computing the exact distribution whenever 

p and 0 , 0 are replaced by their consistent estimators, 
x y 

In spite of this shortcoming, equation (4.5) is used. 

However, Stephens (1976) showed that the asymptotic 

approximation given by Anderson and Darling is conserve- , 

tive as compared to his fitted distribution in the family ; 

of normal distributions [Stephens (1976) Table 4, p. 367] 
and the extreme value distributions [Stephens (1977) j 

Table 1, p. 687]. Thus, it seems reasonable that I 

equation (4.5) is also conservative, that is, if <* n (p) j 

4 

is the true value for the l.h.s. of equation (4 5), then I 

j 

a(p) < a n (p)* ! 


Table 2. Type 1 Errors for Unspecified 


Y 

n 

p = . 25 

.50 

. 75 

1 

20 

1.00 

1.00 

.65 


50 

1.00 

.87 

.25 


100 

1.00 

.50 

.06 

2 

20 

1.00 

1.00 

.55 


50 

1.00 

.74 

. 18 


100 

1.00 

.41 

.04 

3 

20 

1.00 

1.00 

.53 


50 

1.00 

. 70 

. 16 


100 

1.00 

.37 

.03 


29 - ~ / 


From Table 2, it follows that the procedure is only 
sensitive to p whenever p =» . 75 and n > 50. This observa- 
tion was also supported in the simulation study. 

In order to determine the sensitivity of the test to 
the scale parameters, the distance given by 

V s ) " n /o /2 {G(0) " c < e > )2 dG < e > (4.6) 

where G(o) is given in (3.2) and C(o) is given by (3.2) 

whenever tane = sr, s = B /0 = . 90 (. 02) 1 . 10 . Errors in 

x y 

either of the scale parameters can be considered by varying 
s in (3.2). Table 3 contains the expected type I errors 
given by 

Pr[W n > D n (s)j = a n (s) (4.7) 

for different values of n, y, and p. 

From Table 3 one observes that the procedure appears 
to be resilient to errors in the scale parameter and that 
one might have a type I error when y = 3, p = .75, and n = 
100 at the 95% significance level. In addition it also 
appears that the results are symmetric about s = 1. 

5. CONCLUSIONS AND SUMMARY 

A procedure is outlined for determining whether 
a four or five parameter bivariate gamma model is appro- 
priate. The procedure was evaluated and three 
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1 Errors for 
ecified Scales 


CM 

o\ 

. .94 

1.06 

1.08 

1.10 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

.87 

1.0 

1.0 

1.0 

1.0 

1.0 

.96 

1.0 

1.0 

1.0 

.68 

1.0 

1.0 

1.0 

1.0 

.73 

.57 

.89 

.83 

.60 

.40 

1.0 

1.0 

1.0 

1.0 

1.0 

.79 

1.0 

1.0 

.80 

.57 

1.0 

1.0 

1.0 

1.0 

,80 

.64 

1.0 

1.0 

. 66 

.46 

.81 

1.0 

1.0 

.84 

.61 

.46 

.74 

.76 

.49 

.31 

.49 

.78 

.80 

.52 

.34 

.22 

.44 

.46 

.24 

.11 

.98 

1.0 

1.0 

1.0 

1.0 ‘ 

.55 

.86 

.88 

.57 

.38 

.77 

1.0 

1.0 

.80 

.57 

.43 

.70 

.72 

.45 

.28 

.57 

89 

.92 

.60 

.41 

.28 

.52 

.54 

.30 

.16 

.32 

.57 

.59 

.34 

.19 

.10 

.27 

.28 

.12 

.04 


ere omitted since the Type 1 
rs . Likewise whenever n = 20 . 


different functions were evaluated in order to determine 
the procedure's feasibility and sensibility to the nuisance 
parameters. Admittedly, the evaluation is very limited 
and there are several limitations which would prohibit 
this type of procedure as an omnibus test of fit. How- 
ever, the results appear to be promising to the experi- 
mentalist interested in obtaining insight into the stated 
problem. 

There are several nonparametric procedures for test- 
ing (1.2) versus (1.3) and perhaps these are not as sensi- 
tive to the nuisance parameters. However, the proposed 
procedure is based upon "measuring" significant departures 
of the parametric distributions function which are vital 
to the modelers' primary objective. 


RT 2 


- Si J 


} 6 . REFERENCES 

f * 

' I 

\ \ 

[ 1. Anderson, T. W. and Darling, D.A. , (1952). Asymptotic 

| theory of certain goodness-of-f it criteria based upon 

i stochastic processes. Annals of Math. Statist., 23, 

j pp. 198-212. — 

2. Cherian, K. C. , (1941). A bi-variate correlated gamma 

type distribution function. J. Indian Math. Soc. , 5, 
pp. 133-144. " 

3. Darling, D. A., (1956). The Cramer- Smirnov test in 
the parametric case. Annals of Math. Statist . , 26 , 
pp. 1-20. 

4. Gunst, R. F. and Webster, J. T., (1973). Density func- 
tions of the bivariate chi-square distribution. J. 
Statist. Comput. Simul . 2, pp. 275-288. 

5. Jensen, D. R. , (1970). The joint distribution of quad- 
ratic forms and related distributions. Austral. J. of 
Statist . , 12, pp. 13-22. 

6. Johnson, N. L. and Kotz, S., (1972). Distributions in 
Statistics: Continuous Multivariate Distributions' , John 
Wiley, New York. 

Kibble, W. F. , (1941). A two variate gamma type dis- 
tribution. Sankhya , 5, pp. 137-150. 

Mardia, K. V., (1970). Families of Bivariate Distri - 
butions . Hafner Press, New York. 

McAllister, P. R. , (1983). A system for computing 
joint probabilities from Jensen's bivariate F-distri- 
bution. Communications in Statistics, B, 12(2), 223- 
238. ~ 

McKay, A. T. , (1934). Sampling from batches. J.R.S.S. 
B., 1, pp. 207-216. 

Moran, P. A. P. , (1969). Statistical inference with 
bivariate gamma distributions. Bioraetrika, 56, pp. 
627-634. ~ 

Smith, P. E. and Adelfang, S. I., (1981). Gust model 
based on the bivariate gamma distribution. J. Space- 
craft, Vol ■ 18, No. 6 , 545-549. 


1 

A 

it 

r* 

> t 

e 


7. 

8 . 
9. 

10 . 

11. 

12 . 


5* 

M 

1} 


■» V 


4-urtiNiv i\v/i FILMED 


\ 


* 


13. Smith, 0. E., Adelfang, S. I., and Tubbs, J. D. , 
(1983). Computational results for the five-para- 
meter bivariate gamma distribution with application 
to Gust Modeling. Under review Communications in 
Statistics , B. 

14. Stephens, M. A., (1976). Asymptotic results for 
goodness of fit statistics with unknown parameters. 
Annals of Statistics , 4, pp. 357-369. 

15. Tiku, M. L. , (1965). Chi-square approximations for 
goodness of fit statistics U and W . Biometrika , 

52, pp. 530-3. 

16. Tubbs, J. D. and Smith, 0. E. , (1983). A note on 
the ratio of correlated gamma variates. Under review 
by Communications in Statistics , A. 







r 


f ®84 20307 



A DIFFERENTIAL EQUATIONS APPROACH TO THE 
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ABSTRACT 


Analytical and numerical computational methods are 
given for determining the location of the mode as a func- 
tion of the parameters of a class of the bivariate gamma 
distribution. 


I. INTRODUCTION 


Smith, Adelfang, and Tubbs (1983) derived some compu- 
tational results for a family of bivariate distributions. 
In their paper they consider the location of the mode as a 
function of the shape parameters, and \2» and the 
dependence coefficient n . The purpose of this paper is 
to consider this problem in greater detail. That is, the 
paper will consider analytical and numerical computational 
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methods for locating the modal values for the class of density 

functions given in Smith, Adelfang, and Tubbs (1983). The 
general density function is given by: 


'r 1 ’!' 1 -<h +t 2> 

*"1 c 2 ® «> «• 

1 ( t-Ji t-O * Y 1 tf Oi^) " " Z " 1 E 0.. 

(l-,) n r( Yl )r( Y2 - Tl ) J’ 0k “0 Jk 

n^ +k r (Y2-Yx +k ) 


( 1 . 1 ) 


where a. 


jk (l-n) 2 j+k r(Y 2 +j+k)j:k: 


C 1 “ 6 1 X ' c 2 = e 2^' 6 1 ,5 2 are sca l e parameters, Y 2 > Y^ > 1 
are shape parameters, and 0 < n < 1 is associated with the 
correlation coefficient p by the equation n = p * / Y 2 /y ^ • We 
will assume without loss of generality that 53 & 2 ~ 1* 

We will concentrate on the special case = Y^ = Y of 
(1.1) for which the distribution function reduces to 


l(tx* t 2 :Y.n) = 


(t L t 2 ) 


,_X “ (^x’^'^2^ C 1 — n ) 


n ^ ( c xt 2 ) 


j 


(l-n) Y r( Y ) 

This is the form given by Kibble (1941) 


j-0 (l-n)^-’r ( Y +j ) j ' 

( 1 . 2 ) 


Smith and Adelfang (1981) used the above class of density 
functions in modeling wind gust data for the ascent flight of 
the Space Shuttle. A parametric model was selected in that the 
parameters are used to establish engineering constraints for 
the shuttle payload system. Thus, the modal location and value 
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were of interest to this particular application. The authors 
are not aware of any other results, either analytical or numer- 
ical, for the modal location for non-Gaussion multivariate dis- 
tributions. The closest related work is in the area of density 
and mode estimation [3.g. Sager (1978, 1979), de Beauville 
(1978), and Eddy (1980)]. 

In Section 2 we will derive some qualitative results con- 
cerning the behavior of the modal location of (1.2) as a func- 
tion of (y.n). Section 3 presents analogous results for 
another borderline case *» 1, L 2 of (1.1). In Section 
4 we present some numerical procedures based on the theoretical 
investigations of the previous sections. The general case 
Y 2 > Yi > 1 is considered in Section 5. We present some numer- 
ical tabulations for the modal location of (1.1) as a function 
of (y-j,, y 2 > r >) and consider some numerical interpolations from 
the borderline cases considered in Sections 2 and 3. 

2. EQUAL SHAPE PARAMETERS - ANALYTICAL METHODS 

Lemma 1 . The function f(t^,t 2 ;Y.n) defined by (1.2) attains 

O 

its maximum in the region R* = { (t^ , t 2 ) . t^O, t 2 ^ 0 ) on the 
line t-^ = t 2 * 

2 

Proof: Since f is integrable and continuous over R + , it is 

2 

clear that f attains its maximum on R + . Choose any constant 
c > 0. Let h(t) “ f (t,c-t;Y ,n) , 0 < t < c. Then from (1.2) 
we have 
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h(t > - .i 

j=0 J ' • 

where K j (y , n ) > 0 is independent of t. Therefore, 

h ' (t) ■ jf Q K j (*>1) (Y+j-l)C Y+ J" 2 (c-t) T+ -i~ 2 ( c - 2t: ) - 

Since h'(t) > 0 for 0 < t < c/2 and h'(t) < 0 for c/2 < t < c, 
h attains its maximum at t = c/2.. Therefore f^.t^ attains 
its maximum along any line tj+t^c at the point (c/2, c/2). 

This completes the proof. 

Define g(t; Y , n ) * f(t,t; Y ,n). Then by Lemma 1 it is suf- 
ficient to find the point on t > G at which g attains its maxi- 
mum value. Using (1.2) one can show that 

g(t; Y ,n) = c( Y , n ) e ~ 2t/(1 " n) h(t) (2 . 1 ) 

where c( Y , n ) = [ (l~n) (*^n) Y ~ 1 r ( Y ) ] _1 , and h(t) - t^I ^(p^jt), 
where I p (z) denotes the modified Bessel function with* index „ , 
and p(n) = 2/VT7(l-n). 

Using [Abramowitz and Stegun (1964), Eqn. 9-6-28] it is not 
difficult to show that h*(t) » p(n)t Y_1 I y _ 2 (p(n)t), therefore 
f(T,T; Y ,n) is the mode at the bivariate gamma distribution given 
by (1.2) if and only if g’(r) = 2g( T )/(l- n ) or 

*^ I Y - 2 (P( T1 ) T ) = ^-iCp (i)t) , (2.2) 

where p( n ) - 2/^/(l-r,). 

With the aid of (2.1), we may prove the following theorem. 

Iheorem_l. For fixed Y > 1. let t(„> denote the value at which 
f<*(n), * (n) ; Y , n) is a maximum. Then r is continuously differ- 
entiable for 0 < „ <1 and satisfies the initial value problem 


j 

.4 

i 


\i 


I 

- i 

i 


n 

(3 

$■* 
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t ' (n) - (r/2n) ((2t- 2Y+3)" 1 - <l+n)(l-n) l ) 

(2.3) 

T (0) = Y“l- 

Proof. It is easy to show directly from (2.1) that g attains 

its maximum at t = Y - 1 when n = 0, so that t (0) = y “ 1. 

Furthermore is continuously differentiable for 0 < n < 1 and 

computation shows that 3 — 2. / 0 at t = Y”1 and o=0. Therefore, 

at 2 

t ( n) is continuously differentiable in a neighborhood of n = 0 
by the implicit function theorem. The proof will be completed 
by differentiating both sides of (2.1) with respect to 
After some simplification and solving for t * (n ) this yields 

t'(o) = (l-n)I ? (p (n)f )/(4rjq(n) ) - (1+n ) t/2n (1-n ) (2.4) 

where q ( n ) = I^_i(p(n)x) - fr\ I Y _ 2 ^ ^ T ^ * 

By [Abramowitz and Stegun (1964) , Eqn. 9-2-26] , 

I Y _p(p(n)T) = 1^ _ 2 ^P (n ) t ) - (y-1) i y -]_ (p(n ) t )/p(n) t 

l^_l(p(n) T ) = I y _ 1 ( p ( m ) t ) + (y-2) _ 2 (p(n) t )/ p(n)x ? 

substituting these expressions into q(n) and using (2.2) yields 
after some simplification 

q(n) = (1-n) (2 t-2y + 3)I^_ 2 (p(n)T)/2x. 

Substituting this expression into (2.4) completes the proof 
of Theorem 1. 

The nonlinear differential equation (2.3) cannot be solved 
in general in closed form. Some numerical solutions are given 
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in Section 4. However (2.3) does give information regarding 
the qualitative and limiting behavior of t (n) for y > 1. In 
the special case y ~ 3/2, (2.3) reduces to a linear differcn- 
tial equation which can be solved directly by standard methods. 

Corollary 1. If y = 3/2, then 


t (n) = In f mn) 

4/ n l l-/n / 


This result can also be obtained directly from (2.2) using the 
fact that 1^ (z) can be expressed in terms of hyperbolic func- 
tions when y = ±1/2. 

Since the differential equation (2.3) is singular at 
n = 0, its numerical solution requires some additional knowl- 
edge of the behavior of f(n) near n = 0. This is provided 
by the following corollary. 

C orollary 2 . The function t(h) is continuously differentiable 
at n = 0 and satisfies 


t’(0) = -(y~1)/y . Y > 1. 


(2.5) 


Proof . The continuous differentiability of - at n ~ 0 was 
considered in the proof of Theorem 1. Choose n > 0, then by 
the mean value theorem there is a number £ c (0,n) such that 
T (n) = t( 0) + nt'(C) = y " 1 + ri t ' ( £ ) . 

Substituting this expression into (2.3) and simplifying yields 

T ' ( n ) = (y-1+ot ' (Q) (l-f(Hn ) T ' (Q) 

(l-n)(l+2 n T '(£)) 
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Letting n -► 0 and using the continuity of i'(n) wu have 

t ' ( 0) = - (y-1) (1 +t ' (0) ) . 

Solving this equation for x'(o) yields (2.5) 

We will write t(ti,y) when we wish to emphasize the uepend- 
ence of the modal location on y. Theorem 1 and Corollary 2 
may be used to obtain several of the qualitative and asymptotic 
properties of the function t (n/Y) in the region 0 < n < 1/ 

Y > 1. These are summarized m the following theorem. 


Theorem 2, The modal location function x(n#r) has the following 
properties : 

(i) x(n,y) is a decreasing function of n for fixed y > 1* 

3 

(n) lim t (n / Y ) = max (y-^.O) for y > 1; 
n-*-l z 

. . 3 

(m) t ( n / y ) - (y-j) is a decreasing function of y for 
fixed n e (0,1) and y > 1; 


(iv) lim t (n » y ) - (y-~) = 


2 ' 2(l+n) 


, for 0 < n < 1. 


Proo f : We will show that t' (n) < 0 for 0 < n < L Suppose 
not, then since t' ( 0) < 0 by Corollary 2, there is a point 
£ > 0 such that r ' (£ ) = 0 and t * ( n ) < 0 for 0 < n < £. Let 


w(n) = t (n) - (y~ j) and z(n) = 
easy to see that 


- i~n 


2 (1+n ) 


•, then from (2.3) it is 


i , < _ t (n) r 1 1 1 

T (n) 4n w (n ) z (n ) f 


so x'(£) = 0 if and only if w(£) = z{£). 
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Let h = w - z. Note that z'(n) 13 -(1+n) ^ so that h'(0) ** 

w' ( 0 ) - z’ (0 ) - t'< 0 ) + 1 > 0 and h( 0 ) - w( 0 ) - z(0 ) - 0 . 

Therefore since h(e) =■ 0 and h(n) > 0 for 0 < n < S# we must 

have h'(c) 4 0. However, h'(e) = w' (O - z % (c) = t r (O - 

z'(0 = -z' (£) = (1+0*^ > 0. This contradiction proves (i) . 
Furthermore, we have that w(n) > z(n) for 0 < n < 1. 

We will now consider the proof of (iii) . Fix > y 2 > ^ 
and let f(n) » w(n,Y]_) ~ w(n,Y 2 ) where as before w(n,v) “ 
t(ti,y) - (y - ^)- We wish to show that f(n) < 0 for 0 < n < 1 . 

Clearly f(0) - 0 and by (2 . 5) f ' (0) » i - £ <0. Assume 

Y 1 y 2 

to obtain a contradiction that there is a point Ce( 0 ,l) such 
that f(0 = 0. If, in addition, we assume £ is the first 
such point, then f(n) < 0 for 0 < n < 5 so f'(C) >. 0. However, 

using ( 2 . 6 ) at both y^ and y 2 and tbe fact that w(C,y^) “ 

w(£ , y 2 ) it is not difficult to show that 

(yi-y 2 ) 


f’CO - 


4n 


1 _ 1 
w( 0 ~ “ z(;) 


Since y ^ > Y 2 and w(r,) > z(c), it follows that f ' (?) < 0. 
This contradiction completes the proof of (iii) . 

Now we turn to the proof of (ii). First consider the 
case 1 < y 4 3/2. Since t is decreasing in n and positive 
for 0 < n < 1 we know that t* = lim -r(n) exists, where the 
limits at 1 are always from the left. Assume to obtain a 
contradiction that t* > 0. Then it is not difficult to show 
using (2.3) that 
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T ’ (n) 

Therefore, for j < n < 1 we have 

Integrating both sides of this inequality from j to n yields 
T (n) < t C^) + j In(l-n) 

for 1 « n < 1 * However, this implies that T ( n ) - - » aS 
1 + 1. a contradiction. 

The case y > | follows easily from (iii) and the proof of 
(i) because for y ^ ~ 

z(n) < T(n) - (y-|) < t(d,|) 

and both z(n) and t(h,|) approach zero as n •> 1. 

Finally, we consider the proof of (iv) . Let u( n , Y ) = 
w(n,y) - z(n) for 0 < n < 1 and Y > . From the proof of (i) 

we know that u(n, Y ) > 0. From (2.o) we obtain 

w’(n,Y) = 2 ~ fl - + ( 112L f 2(n) - wCn.vl 1 

n L 2 (nj J (_ w (n , y) z Cn^ J 

so that 

w ’(n.Y) < - (Y-|)u(n, Y ). 


u'(n,Y) - w ' (n ,y ) - z’( n ) < - (y -|) u(n , y) + 1. 
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From this inequality we obtain 

^ [u(,. y )e<*>] < e ( ^2 )n 
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u(r, fY )e (Y -2 )rl - u(0, Y ) < — L- I e^ Y '2^ n -1 

(y-^) l -* 

Therefore 

0 < u(n.v) < — K-. 

(Y-|) 

This implies that u(h,y) **■ 0 as Y *► °> and completes the 
proof of Theorem 2. 

3. UNEQUAL SHAPE PARAMETERS- -THE CASE Y L = 1 

In this section we consider another "borderline" case of 
the general bivariate gamma distribution, the case Y] _ = 1. 

For technical reasons we will limit our discussion to the 
range y 2 > 2 and for brevity let y 2 = Y • Then the function 
given by (1.1) reduces to 


f(t, ,t 9 ;l,y,n) = 


t 0 Y -V s 2 


z e“ s l C. 


' j JT 


(3.1) 


where 


(3.2) 


c. - z s j+k r(Y+k- l) . n . , 

J k=0 ^ klr ( Y +j+kT ’ J 0*1.2, .... 

C 1 

and where S]> = ^ s 2 - T ^, and S3 = n s 2 . 


The following lemma allows us to restrict our attention 
to the line t-^ = 0. 


0. 



\ ■ 
\ 


I 


4# 

r ~i 


t 1 


! i 


i- ! 
' i 

f ' ! 
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Lemma_2. The function f (t p t 2 ; l.y, n ) given by (3.1) for 
Y 2, takes on its maximum value in the region t-^ > 0, 
c 2 £= 0 on the line t-^ = 0. 

Proof: Since f is continuous and integrable in the first 

quadrant, we know it takes on its maximum value at some point 
U i* t 2‘ ) - We wil1 P rove that - 0 by showing that for any 
fixed t 2 > 0 , f(t 1 ,t 2 ) is a decreasing function of t x . This 
is equivalent to showing that the function 

S(s) = z e“ s Cj ^ 


j=0 


“j JT 


is a decreasing function on s > 0 where Cj is given by (3.2) 
Note that 


A 


g’(s) = -e" 


ii j_1 


2 yj- + e s i c. ■1^-r 

j=Q J Jl j=0 J jf 


- - e_s ^ c i fr + e ‘ s = c -+i 4 
j ji J=0 j+iji 


j-o 


- -- i 4 (c j" c j+i ) - 

j=0 


Therefore g'(s) < 0 for s > 0 if c j+1 < c for j = 0,1,2 
To this end note that 

c “ „j+l+k r( Y -fk-l) 

k =0 3 kl r (y+j+l+k) 


= r „j + k r(v+k-2) 
k=l 3 Tk'-"l)'!r (y+J+k)’ 


= Z s 3+k r (y+k-1) # k 

k=l 3 klr (Y+J+k)" Y+k^T c j 


l* 
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. k 

since y > 2 implies that 1 1 for k = 1,2,3,.... This 

completes the proof of Lemma 2. 

According to the preceding lemma, the mode of the bivariate 
gamma distribution in this case is the point (0,n) where n is 
the point on t >, 0 where the following function is a maximum: 


g(t) = t Y_1 h(t) 


where 


u/-\ _ T /-nt N r(v+k-l) 
h(t) k ^ 0 ( I=n> k'l r ( Y '+k)' 


z (r%) 


k=0 


1-n kl (y+k-i) 


Note that 


dt ( 


rUL./' 1 


((ffr) h (t) ) = 


j " r Y+k-1 i 

3E ( k y*-> kl(Y+k-l) ) 


. n : ,.t / +k - 2 1 

- r^r k f 0 W ET 


_ n 


/nt_\ Y 2 nt/(l-n) 
-n ''T r n v 


Therefore, 

h(t) - ^ /* ^ y ' z 0 ns/a-„) ds 


s Y “ 2 Cl-n ) ds 


so that the function we wish to maximize is 


g(t) = e -t/(1 - n) I l s Y " 2 e ns/( l‘ n) ds ' 


(3.3) 
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— - " Let y(n) ' or when necessary p( n . Y ). denote the value 
for which f (0,p (n ,y) , 1, Y , ,j) is a maximum where f is defined 
by (3.1) and (3.2). Then p (0) - Y -l, „ („) > Y _ 2 for 0 < „ < 1. 
and satisfies the equation 


S(v) - (1-n )u y ~ 2 e~ u , 0 < n < 1, 


(3.4) 


where g is defined by (3.3). 

Proof: It is easy to see from (3.3) that g attains its maximum 
on [0,-) at a point t* > 0 for which g’ (t*) - 0 and g"(t*) < 0. 
Differentiating (3.3) we obtain 


g,(t) = " I^T 8(t) + t Y “V c 


and 


g ' ,(t > * - l3=n «*(*=> + ^■V t ( Y -2.t), 

Therefore g*(p) - 0 implies (3.4) and g"(p) < 0 implies that 
^ > Y - 2. Since when „ - 0, g(t) - e’ 6 ^ it is easy to 
see that p (0) . Y -l. This completes the proof of the lemma. 

With the aid of these preliminaries we may prove the 
following theorem in the spirit of Theorem 1. 

Theorem 3. For fixed y > 2, let p(n) denote the value of 
which f (0, p (n) , l, Yt n) is a maximum. Then p is continuously 
differentiable on 0 < rj < 1, p ' ( 0 ) = -1 + A, and on 0 < n < 1 
p satisfies the initial value problem 



v'(n) 


_i‘ (»* " (y-1) + nl 
n 0.-n)(u-(Y--i)) 
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(3.5) 


Proof : As in the proof of Theorem 1 the continuous differen- 

tiability of p in a neighborhood of n“0 may be proved by applying the 
implicit function theorem to (3.4). This differentiability will 
be extended to all of [0,1) by proving that (3.5) holds. Let 

g(t,n) denote the function defined by (3.3) and let and g 

t n 

denote its partial derivatives with respect to t and n» 
respectively. Then differentiating both sides of (3.4) with 
respect to n we obtain 

gj-Oj.TOu' +g n (p,n) = (l-n)e" y M Y " 3 (y-2-p)p '-e _lJ p Y “^. (3.6) 

By definition g t (u,n) = 0 and direct differentiation of (3.4) 
and integration by parts yields for 0 < n < 1 that 

s 4 (t ' n) = "7TT? 8(t ’ n) 

(1-n) 


( 1-n ) ' 


■t/ (1-n) f t s y- 1 e ns/(l-n) ds 


— ~ — j g(t,n) 

(1-n) 2 


+ 1 e ~ t/(l-n) j Izjl t Y-l e nt/ (1-n) 

(1-n) 2 n 


- ^ (y - 1) si s Y - 2 e n5/(1 -^ds] 


til • a, 
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" "(l-n) 2 8<t,,1> + tV ' le ' 1 - 8(s.n). 

Therefore, using (3.4) „ e obtain 


E n ( "’ n> ' -1^ i,Y ' 2c ''' + ^^7 » Y 'V" 

- » Y ‘V* 

“ - T '*«-<-Ar + iT7S=5T - 

- uT-2 e -u u-(y-l) 
n 

Substituting this expression into (3.6) and simplifying yields 
(3.5). For n = 0, an easy calculation shows that 


£ n <t: ' 0) = ~7xhj 


e _t t Y 


from which substitution into (3.6) with „ - 0 and u = Y - i 

shows that u ' (0) = -1 + I TK „- „ „ , 

v J 1 + Y ' Th completes the proof of 

Theorem 3. 

The following corollary exploits the fact that (3.5) 
reduces to a linear differential equation when Y « 2. 

C°l£±±ai2_ 3. If Y - 2, then 


v (n) 


ln ^T^» 0 < n < 1. 


Proof. For Y = 2 equation (3.5) reduces to 

u ' + U _ 1 
n (l-n) n 
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which is easily solved in closed form by standard methods to 
show the desired result. This result is also easily derived 
directly from (3.4). 


It is interesting to note that the translated modal 
location function v( n ) - y( n ) - (y-2) satisfies the differ- 
ential equation 


v‘ <n) 


v(n) + y-2 
n 



v(0) = 1, v* (0) - -1 + i 

Y 


whereas the translated modal location function w(n) = 

* ^ 

‘ °f Section 2 satisfies the analogous differential 
equation 


w’(n) 


w(n) + y - | 
4n 



2(l+n) , 

1-n J 


w(0) - i . w' (0) » -1 + I. 

Z Y 

For this reason ji behaves in a manner similar to t. Its 
properties are stated in the following theorem. Since the 
proof of this theorem is entirely analogous to the proof of 
Theorem 2, it is omitted. 

Theorem_4. The modal location function uCn.y) has the follow- 
ing properties: 

(i) m(ti,y) is a decreasing function of n for fixed 
Y > 2; 

(ii) limp( n ,Y) = y - 2 for T > 2; 
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(iii) p(n,v) - (y~2) is a decreasing function of y for 
V > 2 and fixed n e (0,1); 

(iv) lim(ij(n,Y) - (y~2)) = 1-n for 0 < o < 1, 

Y->~ / 

4. NUMERICAL RESULTS 

( In this section we present some quantitative results based ; 

! on the results of the previous sections. Table 1 shows the j 

! 

j value of the modal location function for equal shape parameters j 

| t ( n» y) for various values of n and y. Table 2 shows values of J 

I 3 1 » 

the translated modal location function w(h,y) = t(n,Y) “ (y"^)* \ 

1 ? 

{ This table illustrates the qualitative behavior of this func- ] 

t i 

tion derived in Theorem 2. The limiting values of n = 1 and j 

f 

I 

Y = ® are taken from Theorem 2. | 

The values in Tables 1 and 2 were computed using Theorem ’ 

1. Specifically, a fourth-order Runge-Kutta algorithm was used 
to compute an approximate solution of the differential equation 
(2.3) on the interval 0 < n < 1 for each specified value of y. 

Since equation (2.3) is singular at ii = 0, Corollary 2 was used 
to replace the initial condition x(0) “ hy the approximate 
initial condition 

/ 

> 

T(h ) = Y -1 - (III) h j 

r 

where h is the step size of the numerical method. Figure 1 l t ‘ 

shows the data of Table 2 in graphical form and illustrates the L * 

f * 

t 

i 

behavior of the function w(d,y) derived in Theorem 2. 




/ 


Tables 3 and 4 show the corresponding results for the 
modal location function w(n,y) for the case = 1, Y£ “ Y 
and its associated translate v(n,y) ■* n(n,y) ~ (y-2). These 
tables were computed by the same methods as Tables 1 and 2 
except using the results of Section 3. Figure 2 illustrates 
the qualitative behavior of the function v(n,y) as indicated 
by Theorem. 4. 

Note that the differential equations (2.3) and (3.5) 
allow the modal location to be computed recursively in n for 
a fixed value of y as a dynamic process in a time scale meas- 
ured by the modified correlation coefficient r, . Error in the 
computation is introduced through the discretization of this 
continuous evolutionary process. A more conventional compu- 
tation of the modal location would require an independent 
calculation for each value of n with error introduced through 
the truncation of the series representation (1.2) of the dis- 
tribution function. This error becomes particularly trouble- 
some as i) + 1, 
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Table 1. Selected values of the modal location function t(h,y) 
for equal shape parameters. 


i 

i 


; 

0 

.1 

.3 

.5 

.7 

.9 

1 

« 

; i.i 

.1000 

.0908 

.0720 

.0525 

.0322 

.0110 

.0000 

: i.3 

.3000 

.2765 

.2268 

.1724 

.1117 

.0412 

.0000 

1.5 

.5000 

.4660 

.3931 

.3116 

.2169 

.0958 

.0000 

2.0 

1.0000 

.9491 

.8412 

.7277 

.6127 

.5279 

.5000 

3.0 

2.0000 

1.9334 

1.8034 

1.6862 

1.5^39 

1.5268 

1.5000 

10.0 

9.0000 

8.9152 

8.7754 

8.6697 

8.5891 

8.5624 

8.5000 


c_n 

OJ 


Table 2. Selected values of the translated modal location function 
w(n,y) for equal shape parameters. 



0 

. 1 

.3 

. 5 

.7 

.9 

1 

i. i 

.5000 

.4908 

.4720 

.4525 

.4322 

.4110 

.4000 

1.3 

.5000 

.4765 

.4268 

.3724 

.3117 

.2412 

.2000 

1.5 

.5000 

.4660 

.3931 

.3116 

.2169 

.0953 

.0000 

2.0 

.5000 

.4491 

.3412 

.2277 

.1127 

.0279 

.0000 

3.0 

.5000 

.4334 

. 3034 

.1862 

.0939 

.0268 

.0000 

10.0 

.5000 

.4152 

.2754 

.1697 

.0891 

.0264 

.0000 

OP 

.5000 

.4091 

.2692 

.1667 

.0882 

.0263 

.0000 


O O 
-n x, 

x Z\ 

O 

S? 

iO '• 
C. V* 

r* 

t_ to 
-1 «.-> 





Table 3. Selected values of the modal location function w(n.Y) 
for y l = 1, y 2 = V > 2. 


y\i 

0 

.1 

.3 

.5 

.7 

.9 

1 

2.0 

1.0000 

.9482 

.8322 

.6931 

.5160 

.2558 

.0000 

3.0 

2.0000 

1.9310 

1.7767 

1.5936 

1.3702 

1.1111 

1.0000 

10.0 

9.0000 

8.9085 

8.7170 

8.5155 

8.3081 

8.1011 

8.0000 


Table 4. Selected values of the translated modal location function 
g v(n ,y) for y ]_ = 1 . y 2 = Y i 2. 


y\h 

0 

. 1 

.3 

.5 

.7 

.9 

1 

2.0 

1.0000 

.9482 

.8322 

.6931 

.5160 

.2558 

.0000 

3.0 

1.0000 

.9310 

.7767 

.5936 

.3702 

.1111 

.0000 

10.0 

1.0000 

.9085 

.7170 

.5155 

.3081 

.1011 

.0000 

CO 

1.0000 

.9000 

.7000 

.5000 

.3000 

.1000 

.0000 
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UNEQUAL SHAPE PARAMETERS 


In this section we briefly consider the mode of the 
general bivariate gamma distribution given by (1.1). By 
setting the partial derivatives of f (t-^, t 2'" , 'l' Y 2' n ^ w *- tk 
respect to t-^ and t 2 equal to zero, one finds that f attains 
its maximum at the point (tpt£) whose coordinates satisfy 


1-n 


CO 00 


ti “ V- 2 2 a ik (Y i + j" 1) 

1 15 k=0 j=0 JK x 


(5.1) 


and 


cn co 


t 9 = 2 Z a, 1 .(Y2 + j +k - 1 ) 

1 b k=0 j=0 z 


(5.2) 


CO 00 


where S = Z Z a., 
k=0 j=0 JK 

and a^ given as in (1.1) depends on t^ and t£. 

Table 5 shows selected values of the modal location for 
the case y£ = 3. They were computed by truncating each of 
the series in (5.1) and (5.2) to about fifty terms and simul- 
taneously iterating on these equations until an approximate 
solution is obtained. These computations become unreliable 
as n + 1 and the truncation error becomes unacceptable. 

Figure 3 gives a graphical representation of the change 
in modal location with n and y ^ for fixed Y 2 = 3. It is 
interesting to note for a fixed n the extent to which the 
modal location may be approximated by linear interpolation 
between the borderline cases discussed in Sections 2 and 3. 
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More specifically, we have the empirical approximations 



T(n.Y 2 ) 


(5.3) 


and 

Y i~ 1 

t 2 - u(n,Y 2 ) + (x (n ,y 2 ) - v(n,Y 2 )) (5. A) 

where x and y are as defined in Sections 2 and 3, respectively. 

This empirical relationship is a subject for further investiga- 
tion. 
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Table 5. 

• o 

f\ 

Location of the mode using (5.1) and (5. 2) with y? ** 3. 

Approximate values computed using (5.3) and (5.4J are 
denoted by *. ■ When these values are equal to two 
decimal places, only one is given. 

o 

\ Y 1 

1 

1.5 

2 

2.5 

3 , 

0 

(0,2.00) 

(.50,2.00) 

(1.00,2.00) 

(1.50,2.00) 

(2.00,2.00) 

.25 

(0,1.82) 

(.46,1.82) 

( .92,1.83) 

(1.38,1.83) 

(1.84,1.84) ? 

.50 

(0,1.59) 

(.42,1.62) 

( .84,1.64) 

(1.26,1.66) 

o 

(1.69,1.69) I' 

.75 

(0,1.31) 

(.40,1.38) 

(.39,1.37)* 

✓-N/-S 

• • 
00 
o o 

1 f— 1 
• • 

4> 

4> 

v-'vy 

X- 

(1.19.1.51) 

(1.18.1.51) * 

i 

(1.58,1.58) 

.85 

(0,1.18) 

(.41,1.26) 

(.39,1.27)* 

( .80,1.36) 

( .77,1.36)* 

(1.17.1.45) 

(1.16.1.45) * 

(1.54,1.54) 

( < 1 . 

> 





i 

' f 

- ! 

1 ‘ 1 
l 1 

1 1 

» C' ^ 

" V j 

4 1 
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ANALYSIS OF WIND GUST DATA 


J. D. Tubbs 

Department of Mathematical Sciences 
University of Arkansas 
Fayetteville, Arkansas 


ABSTRACT 


This paper summarizes the analysis of wind gust data 
using statistical and mathematical procedures which were 
developed for the bivariate gamma distribution. 

1 . INTRODUCTION 

Adelfang and Smith (1981) discuss the use of the gamma 
distribution in modeling gust data at Cape Canaveral, 
Florida. Smith and Adelfang (1981) treated gust amplitude 
and length scale as the variables of the bivariate gamma 
distribution. Smith, Adelfang, and Tubbs (1983) presented 
some useful analytical and computational results for a class 
of the bivariate gamma and applied some of these results to 
the wind gust data. The purpose of this paper is to analyze 
the wind gust data using some additional analytical results 
obtained for the bivariate gamma distribution. 

2 . DATA 

The data used in this paper consists of absolute gust 




62 


magnitude and gust length for both the zonal and meridional 
components. The 150 wind profiles were filtered using the 
band pass filter for wavelengths within 420-2470 meter band. 
Data were available for the reference altitudes: 4Km, 6Km, 
8Km, lOKm, 12Km, and 14Km. The data was paired into bivar- 
iate components for both the zonal and meridional components, 
denoted by the pairs (Au,Lu) and Av.Lv) , respectively. 

3. ANALYTICAL PROCEDURES 

The data were partitioned according to reference alti- 
tudes, then the 150 observations were analyzed using both 
univariate and multivariate techniques. Simple descriptive 
univariate techniques were generated using PROC UNIVARIATE 
in SAS. These procedures were used to help in the assess- 
ment of the marginal distribution. The multivariate 
descriptive procedures consisted of bivariate scatter plots 
and contour plots. 

Goodness of fit tests consisted of a univariate test 
for marginal normality generated by SAS, two tests for 
bivariate normality as discussed in Meredith and Tubbs 
(1981), and a bivariate test for the gamma distribution. 

The latter procedure is a bivariate Chi-square type test 
which uses the computational results for the distribution 
function as presented in Smith, Adelfang, and Tubbs (1983). 

Parameter estimates for the bivariate gamma distribu- 
tion were evaluated. These estimates were then used in 


63 


• a* ^ • • 


I 


generating the three-dimensional bivariate gamma density 
function plots and the modal locations were estimated 
using the results given by Brewer, Tubbs, and Smith (1983). 

4. RESULTS 

The results are summarized in Tables 1-7. Additional 
results are given in Appendices A and B. 

Tables 1-6 summarize the results for both the test of 
fit and the parameter estimates for the bivariate gamma 
distribution. There are two main tests for bivariate nor- 
mality and both of these are discussed in Meredith and 
Tubbs (1981). The first is a procedure proposed by Rincon- 
Gallardo et al. (1979). Since this procedure transforms 
the data to a univariate test for uniformity three differ- 
ent tests for uniformity are used. The second test for 
bivariate normality is based upon a procedure proposed by 
Cox and Small (1978). 

The bivariate test for the gamma distribution is a Chi- 
square type test of fit. Thus, this procedure has the 
usual difficulties of selecting the number of cells and 
cell location that are associated with this type of test. 

In the interest of time and space a fixed procedure was 
applied for all the data sets. Namely, the marginal distri- 
butions were partitioned according to the .05, .25, .50, .75, 
and .90 quantiles based upon the gamma parameter estimates. 
This particular choice affected the results for some of the 
data sets, however, it seemed a reasonable global choice. 
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The univariate gamma parameters were estimated using 
a maximum likelihood procedure presented by Greenwood and 
Durand (1960) and discussed in Tubbs and Brewer (1981). 

Appendix A contains the results for the univariate 
descriptive statistics. Appendix B contains plots for each 
data set. The density functions were generated using the 
gamma parameter estimates. The contour plots are level 
slices of the density function and are not equal proba- 
bility contours. The location of the mode is denoted by 
the symbol + and this value is computed using the analytical 
results given in Brewer et al. (1983). Table 7 summarizes 
the results for the modal location. 
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Table 1. Summary for Wind Gust Statistic 

Using Band Filter 420-2470 Altitude - 4 Km. 


Multivariate Test 


Gamma 


It 

(Au.Lu) 

(Av.Lv) 

Cramer-Von Mises 

.2062 

.2144 

Watson's U^ 

.2023* 

. 2058* 

K - S 

.0618 

.0551 

Cox 

11.45** 

26.07*** 

Chi-Square 

33.3 

53.00*** 


Univariate Test 


Au 

Normality 

Lv 


.077* 
.068 
. 090** 
.077* 


< 

i 

\ 

i 

$ 

i 

< 


Parameter Estimates 


Au 

3.402 

2.430 

.2280 

Lu 

5.275 

.007 

Av 

2.808 

1.891 

.3415 

Lv 

4.429 

.006 


* denotes that test is significant at .05 level. 

** denotes that test is significant at .01 level. 

*** denotes that test is significant at .001 level 
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Table 2. 

Summary for Wind Gust 
Using Band Filter 420- 

Statistic 
•2470 Altitude 

= 6 Km. 

/ 

/ 

t 

Multivariate Test 

(AuLu) 

(Av.Lv) 

Normality 

Cramer-Von Mises 
Watson's U^ 

K - S 

.3806 

.2411* 

.0897 

.2942 

.2208* 

.0623 

t 


Cox 

24.4*** 

31. 8*** 

- 

Gamma 

Chi-Square 

72.03*** 

57 . j.4*** 

t 


Univariate Test 


Au 

Normality ^ 
Lv 


.047 

.081* 

.054 

.072 


Parameter Estimates Y P p 


Au 

2.577 

1.916 

.4217 

Lu 

4.102 

.005 


Av 

3.168 

2.030 

.2506 

Lv 

4.895 

.005 



* denotes that test is significant at .05 level. j 

** denotes that test is significant at .01 level. 1 

*** denotes that test is significant at .001 level. 

I 

\ 

t 
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Table 3. Summary for Wind Gust Statistic j 

Using Band Filter 420-2470 Altitude = 8 Km. \ 

t 

> 

i 


Multivariate lest 


(Au.Lu) 

(£,y,Ly) 

Cramer-Von 

Mises 

1. 090*** 

. 490* 

Normality Watson's 


. 721*** 

. 409*' 

K - S 


. 102* 

.083 

Cox 


10.67** 

8.67* 

Gamma Chi-Square 


34.74 

52 . 57** 

Univariate Test 




Au 


.114** 


Normality ^ 


. 108** 
. 107** 


Lv 




Parameter Estimates 

a 

Y 

A 

s 

A 

P 

Au 

3.023 

2.113 

.3396 

Lu 

4.149 

.005 


Av 

2.922 

1.811 

.4253 

Lv 

4.614 

.005 



/i 



/ 

/ 

/. 


i 




* denotes that test is significant at .05 level. 

** denotes that test is significant at .01 level. 

*** denotes that test is significant at .001 level. 
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Table 4. Summary for Wind Gust Statistic 

Using Band Filter 420-2470 Altitude = 10 Km. 


Multivariate Test 

(Au.Lu) 

(Av.I.v) 

Normality 

Cramer- Von Mises 
Watson's 
K - S 

.129 

.126 

.052 

. 753** 
.469** 
.104* 


Cox 

8./3* 

8.58* 

Gamma 

Chi-Square 

45 . 04** 

36.29 

Univariate Test 




Normality 

Au 

Lu 

Av 

Lv 

.073* 
.063 
. 073* 

. 109** 


Parameter Estimates Y 

A 

b 

p 


Au 

3.041 

1.909 

Lu 

5.203 

.006 

Av 

2.522 

1.300 

Lv 

4.543 

.005 


* denotes 
** denotes 
*** denotes 


that test is significant at 
that test is significant at 
that test is significant at 


.05 level. 
.01 level. 
.001 level. 


i 1 



I 


69 



/ 


t 


Table 5. Summary for Wind Gust Statistic 

Using Band Filter 420-2470 Altitude = 12 Km. 


Multivariate Test (Au.Lu) (Av,Lv) 



Cramer- Von Mises 

. 465* 

.398 

Normality 

Watson's U^ 

.391** 

. 329** 

K - 

S 

.077 

.076 


Cox 


25.04*** 

23 . 70*** 

Gamma 

Chi- 

■Square 

52.54** 

41.24* 

Univariate Test 






Au 


. 116** 


Normality 

Lu 

Av 


. 112** 
. 042** 



Lv 


.085** 


Parameter Estimates 

Y 

l 

p 


Au 

2.155 

.969 

.2983 

Lu 

4.113 

.005 


Av 

2.612 

1.051 

.3066 

Lv 

4.023 

.005 



* denotes that test is significant at .05 level. 

** denotes that test is significant at .01 level. 

*** denotes that test is significant at .001 level 
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Table 6. Summary for Wind Gust Statistic 

Using Band Filter 420-2470 Altitude - 14 Km 


Multivariate Test 



Cramer-Von 

Mises 

1.359*** 

.616* 

Normality 

Watson's U^ 


. . 647*** 

.357** 


K - S 


. 134** 

.084 


Cox 


14.6*** 

17 . 9*** 

Gamma 

Chi-Square 

48 

48.29** 

42.97* 

Univariate Test 






Au 


. 090** 


Normality 

Lu 

Av 


. 096** 
. 085** 



Lv 


.061 


Parameter Estimates 

Y 

& 

P 


Au 

3.803 

1.549 

.2420 


Lu 

4.725 

.005 



Av 

3.324 

1. 161 

.3171 

* 

Lv . 

5.057 

.006 



* denotes that test is significant at .05 level. 

ww denotes that test is significant at .01 level 

denotes that test is significant at .001 level 
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Table 7. Modal Location 


Altitude 

Variables 

4000 

(Au.Lu) 

(Av.Lv) 

6000 

(Au.Lu) 

(Av.Lv) 

8000 

(Au , Lu) 
(Av.Lv) 

10000 

(Au.Lu) 

(Av.Lv) 

12000 

(Au.Lu) 

(Av.Lv) 

14000 

(Au.Lu) 

(Av.Lv) 


Method* 



I 


II 

.939 

.874 

581.9 

516.7 

.939 

. 873 

581.7 

516.7 

.747 
1. 008 

539.6 
129. i* 

.745 

1.007 

539.4 

729.4 

.875 

.965 

571.0 

639.8 

.874 

.960 

570.8 

639.8 

.984 

1.064 

624.8 

624.8 

.977 

1.058 

625.2 

624.8 

1.083 

1.402 

556.8 

546.2 

1.080 

1.400 

556.8 

546.2 

1.713 

1.868 

705.4 
626 . 0 

1.713 
1. 865 

705.4 
626. 0 


Method I 
Method II 


Truncation of a double 
Interpolation. 


series. 
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5 . SUMMARY 

The data sets are discussed according to reference 

altitude. 

4 Km. The normal or the gamma are not rejected for the 
zonal (u) components. As discussed in Meredith and 
Tubbs (1981) the Cox and Small procedure is sensitive 
to symmetry and is not recommended for this data. The 
gamma distribution was rejected for the v-component and 
normality was accepted. However, marginal normality 
was rejected at the .01 level for the absolute gust 
magnitude (Av) . 

6 Km. The bivariate gamma was rejected in both wind com- 
ponents and normality was not rejected. 

8 Km. Normality was rejected for both wind components. The 
bxvariate gamma was accepted in the u-component but not 
for the v-component. 

10 Km. The u-component appears to be normal whereas the 
gamma is accepted in the v-component. 

12 Km. Both distributions appear to be suspect for the 

u-component and the gamma is accepted for v. Normality 
is also rejected for v by considering the marginal dis- 
tributions. 

14 Km. Neither distribution is acceptable for u and the 
gamma is perhaps better for v. 
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6. CONCLUSIONS 


The wind gusc data was analyzed using some new proce- 
dures for the bivariate gamma distribution. The analysis 
was meant to be informative, in that it represents examples 
for some of the analytical procedure. The analysis is not 
meant to be completely thorough. Hence, there are still 
some unresolved questions concerning the applicability of 
the bivariate gamma for modeling wind gust data. One sus- 
pects that neither the normal nor the gamma are completely 
appropriate, however, perhaps both could provide acceptable 
results for defining engineering constraints. 

As mentioned in the paper the test for gamma is a Chi- 
square type procedure which has inherent problems which 
does not lend itself to easy data independent analysis. 
Instead it requires judicious selection of parameters. 

This analysis did not take advantage of this option, hence, 
the rejection of the gamma could be attributable to poor 
cell location choices. 

Every data set was analyzed using a test for equality of 
shape parameters as proposed by Tubbs (1983) . This hypothe- 
sis of equal shape parameters was rejected in every case. 
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appendix a 



H*i variaCe sun ®>ary statistics generated usine 
PROC UNIVARIATE is the Statistical Analysis 
Svstem. ~ — J 




The next four pages contain the summary statistics 
for each of the univariate variables Au, Lu, Av, and Lv, 
respectively. The reference altitude is 


ALTITUDE - 4000 



/ 


y . 




or;r' r ■- n 

OF POOR QUrii-fTY. 


nocnns 


rjy. 

'•'ID 1>I.V 

rirvlts 

y 

:-.i u: -~o 
a r:< ?: 
u. c '^ o 

'riT'i-j j'I, 


3 50 

1.2 err/, 

0.706/16 
0.602612 
360.207 
50.0 C4 
24.2627 
5CG2.5 
ISO 

0.0777153 

D71VSII rs 


LG TEST 
0.05 
0.006 
0 .20 
0.3 
0.33 

tr: leat 

36 7 

30 

32 

30 2336 
:: CO 7 
2C 3 
20 52 

22 00172030556 
20 2LO1357S0 


FUil l r-2£ 

cur: 

vapiw.cc 

::uraor>is 

cos 

i -I-.,. 

prrc’RTc 
rar , 5 vsr 

FOOD’ D 


HIGHEST 

3.12 

2.13 

3.13 
3.16 
3.67 


150 

1001 1 IAX 

3.67 

209.976 

751 C3 

1 .7825 

0.495C23 

501 MU 

1.295 

C.23930C 

251 pi 

0.0775 

74.3505 
. 05767 C 6 

01 mi: 

0.05 

C .0001 

IVJKIE 

2.62 

0.0001 

C3-C1 

0.905 

0 . 02 / 

i:wn 

C.G 8 


ixrPLor 
o 


4 

-> 

u 

1 

2 

11 

( 


3.7+ 

1 


V 


10 12311 

*' u 

r 

V 

16 14777913345089 

* i* 

1 * 

1 .9-t- 

14 12056670S20035G6G802 

20 ' «’ 

rr 

0 

1? 111234C725C77 

13 

• ** 

10 123334567702055C70 

10 < 


5 1334567 L01111’35CC£ 2 

12 ] 


6 23CC8821235629 

14 V 

c 

4 3473G7C99 

2 ' 

% 

2 00356 

•s 1 

5 ' 

Y 

0 52 

u it 

2 *i 

5 


(XJA1ITILES(DEF=4) 


291 

955 

SOI 

101 

51 

11 


3.402ES 

2.863 

2.339 

0.59 

0.3905 

0.0GS36 


IULilPLjf STIJ'.LEAF BY 1C**-01 


0 . 1 t+ * 


1C71IL P.3C2/i'ILlTY PLOT 




** ** 
*vH 
*+ 


-fX*** 

**** 

**** 


***# 


WK ***.,_}. 


*i — — i— - 
-2 


-1 


40 


-+ • 

+1 


'%■»* 


80 



V/ iRIADLE^LU 



m e 

ins 


N 

150 

sun HITS 

150 

LEAK 

600.568 

sun 

121345 

STD DEV 

344.685 

W.M/’Jv'CE 

110008 

SKOJHESS 

0.694665 

injiriosis 

0.642141 

USS 

1150GG718 

css 

17702335 

CV 

42.60C 

STD hear 

20.1434 

T:rEAi:-0 

20.7445 

pros', vrv 

0.0001 

SGI I PARK 

5C62.5 

PROEM, ST 

0.0001 

;u. c = o 

150 



D.-rOEYE 

0.0689676 

rrcr.’.D 

0.000 


QUA] TITLES (DET=4) 


100% MAX 

1974.9 

59% 

1964.6 

75% 03 

1059 

95% 

1415.38 

50% i:ed 

786.45 

50% 

1277.49 

25% Q1 

537.3 

10% 

395.91 

0% I'll; 

88. 3 

5% 

347.845 

radge 

03-01 

1'ODE 

1886.6 

521.7 

702.4- 

1% 

158.625 


EXTRiy.ES 


1950+ 


LOPu/X PPOEABILITY PLCT 


LQ. 1ST 

HIGHEST 



80.3 

1435.1 



226.2 

1480. 7 



‘ 24G.8 

1792.7 



304.7 

1554.7 



334.4 

1974.9 



bra* leaf 


> J ; dgeplo 

19 57 


2 0 


18 




17 9 


1 V 


16 


f 


15 


( 


14 03445 


5 <■ 


13 12345 


5 ; 


12 01789 


5 V 


11 0111333455GG7S 


14 • 5 


10 C01236G6778 


11 ! 

mm 

9 C0112346C9S 


11 T 


8 0112233334556777895 

19 V + 


7 000335566778899 


15 * 

-* 

6 122245CC5959 


12 1 

rr 

5 001223 4445C6677777G55 

21 -i 


4 000122245G777C9 


15 


3 03455567759 


11 V 


2 35 


2 


1 


ft 


C 9 


1 


8 * ' 

“ “» 




IULTIPLY STS . . LETT DY lG^-102 


‘J 

U 

1! 


tt 

H 

V 
u 

? 

V 

M 

T 

<\ 


S 


* 

***** 


50 i * 


•f*** 

***** 

.,-*** 

*** 

**y ** 

***** 

********.$. 

* * - 14 + 

-H -f 

+ n 




-V:r.IAELD=AV 


OF POOR QU/vLtTY 


miEHS 


OJ/JnTLCG(CLr=4) 


II 

150 

cui: ijgts 

150 

100% IIAX 

5.96 

991 

4.93488 

rusAi: 

1.4844 

cun 

222.66 

751 Q3 

1.86 

951 

3.012 

ltd DEV 

0 .822324 

VARI/FCE 

0.G7G217 

SGI I1ED 

1.405 

90% 

2.4G1 

Fn: ’IIECf 

1.59111 

1URT3CIC 

5.G5C1G 

251 01 

0.9 

10% 

0.673 

PCS 

431.273 

css 

100.75G 

0% hill 

0.03 

5% 

0.51 

cv 

55.3077 

era ;■£ m 

0.0671425 

• 


1% 

0.03 

T:.';L7i'=0 

22.1082 

PROD’. IT V 

0.0001 

PVii JGE 

5.93 



lg.: r</:n: 

5GG2.5 

PROL’/„S'! 

0.0001 

G3-C-1 

0.96 



‘ .1" I c = 0 

150 



HCDE 

1 



3:I50R.JAL 

0.0905899 

PRCCl.sD 

1-0.01 






DiraEilES 


LOI.TST 

0.03 

0.03 

0.07 

0.1 

0.15 


HIGHEST 


3.54 

3.54 

3.8 

3.05 

5.06 


'Y.l . LEAF 
G 0 
5 
5 
4 
4 

3 55C9 
3 11 

2 5555G709 
2 00011111222222334 

1 55555555C6666G6G777777777o0CC0f8t , 999599 
- 1 000000000900111111111222222233333444444 
0 553C0GG677777CC0C0C00oSC8099£55999 
; C 001112 

‘ 6.25+ 

ft 

»• 


i r.cr?r,ar 

l 


4 0 

2 ? 

C 

17 5 

39 + + 

39 * .j * 

34 + -t 

6 S 

1JCHHAL PnOEADILITY PLOT 


4 

I 

4 


\ 


j 


< 


t 

f 


* 




J *.* * 
**'!.. J 


+***++ 

******** 

******* 


************ 
0.25+*** ****++++++ 

+ -2 -1 


+0 


”+l H +2 


I 


i.orn.Ts 


C-U/j .TILES (DCT=4) 


i; 

150 

SUI3 \JGTS 

150 

100% MAX 

2017.1 

99% 

1065.52 


IJ713 

723.427 

sun 

10C514 

75°. 03 

932.275 

95% 

1315.47 

\ 

\ / 
\ * 

STD DEV 

329.434 

VARIAl’CE 

108527 

50?; ilED 

C67.1 

905 

1152.19 

SEIMJESS 

0.709547 

wWiOGIS 

1.17774 

25?. 01 

500.975 

10% 

350.09 

l ' 

DOS 

94672440 

CSS 

16170518 

0% HIM 

96.1 

5% 

244.26 

/ 

/ 

cv 

45.538 

STD 3J713 

26.8982 



1% 

104.566 


T:iXAi:=0 

26.895 

PROB’.'SlTV 

0.0001 

RADGE 

1921 




ogii n/m: 

5662.5 

pncn.’/,;S'f 

0.0001 

G3-Q1 

431.3 




im; c =-- o 

150 



1'iCDE 

601 .1 




D si 302371. 

0.0776302 

PR031.D 

0.024 







d:tr!:;ies 


la.’EST 
96. 1 

112.7 

122.7 
134 

167.1 


HIGHEST 

1427.3 

1477.5 

1490.1 
1719.9 

2017.1 


OF POOR QUALITY 


1309 

6 

567780 

45558 

0237779 


i; E033PLOT 

1 0 


1 0 


011344566709 

002222334555709 

00C001123333445667777595 

00222233344445677000 

1122233333444C79S 

12235556CC 

270 

012377 


i ULTIPLY STEH.LEAT BY 10 **+02 


■ 2 1 * — - 


1050-1 

V 


kopjtl pnczvxiLiTY plot 


<*** 

*'i~rH 

***+ 

**-(. 

+*** 

***** 

**** 

+*** 

***** 

***** 


250+ **+++ 

1! * ****++ 

50+* + I " ! - 1 

-1 1 -( H -i -i 

-2 -1 +0 
83 


+1 +2 



The next four pages contain the summary statistics 
for each of the univariate variables Au, Lu, Av, and Lv, 
respectively. The reference altitude is 

ALTITUDE ® 6000 
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e next four pages contain the summary statistics 
or each of the univariate variables Au, Lu, Av. and Lv 
respectively. The reference altitude is 
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The next four pages contain the summary statistics 
for each of the univariate variables Au, Lu, Av, and Lv 
respectively. The reference altitude is 

ALTITUDE - 10,000 
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The next four pages contain the summary statistics 
for each of the univariate variables Au, Lu, Av, and Lv, 
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The next four pages contain the sunnaary statistics 
for each of the univariate variables Au, Lu, Av, and Lv, 
respectively. The reference altitude is 
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